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ABSTRACT 

A  new  algorithm  is  proposed  for  the  reduction  of  the  complexity  of  observed 
seismic  waves  in  the  1-100  sec  period  band.  The  algorithm  transforms  the  waveform 
for  an  event  into  a  simplified  form  and  is  referred  to  as  the  Earth  simplifying  transfor¬ 
mation  (EST).  The  event  is  referred  to  as  the  primary  event,  for  which  properties  of 
the  source  are  to  be  determined.  The  EST  algorithm  uses  data  and  synthetics  for 
another  event  nearby,  which  is  referred  to  as  the  secondary  event  The  synthetics  are 
computed  for  a  given  Earth  model  using  a  priori  source  parameters  for  the  secondary 
event.  Data  for  the  algorithm  are  regional  or  teleseismic  seismograms,  including  sur¬ 
face  waves  and/or  body  waves.  In  this  algorithm,  data  for  the  primary  event  are 
deconvolved  with  the  residual  waveform  which  is  obtained  from  the  deconvolution  of 
data  and  synthetics  for  the  secondary  event.  The  residual  waveform  contains  the 
effects  of  Earth’s  structure  on  seismic  waves  which  are  not  predicted  for  the  Earth 
model.  The  simplified  waveform  obtained  from  the  algorithm  is  referred  to  as  the  EST 
seismogram.  The  EST  seismogram  represents  the  waveform  which  is  expected  for  the 
seismic  waves  from  the  primary  event  for  the  Earth  model.  The  source  properties  for 
the  primary  event  can  be  determined  from  EST  seismograms  instead  of  original  data, 
which  contain  not  only  the  information  on  the  source  but  also  the  effects  of  Earth’s 
structure  unpredictable  for  the  model. 

Several  tests  were  made  using  data  from  earthquakes  located  in  southern  Califor¬ 
nia  recorded  at  Global  Seismographic  Network  (GSN)  stations,  with  the  distance 
between  the  primary  and  secondary  events  being  less  than  100  km.  These  tests  indi¬ 
cate  that  EST  seismograms  provide  superior  resolution  power  in  the  source  characteri¬ 
zation  than  original  seismograms.  The  major  limitation  of  the  procedure  is  that 
a  priori  information  on  the  source  of  the  secondary  event  is  required.  This  procedure 
can  be  applied  to  characterize  an  explosion  source  if  source  parameters  for  an  event 
nearby  are  well  known,  with  the  secondary  event  being  either  an  explosion  or  an  earth¬ 
quake.  The  application  of  this  procedure  may  be  useful  for  improving  the  regional 
monitoring  and  discrimination  capability  under  a  potential  CTBT  as  well  as  basic 
seismological  knowledge  on  source  properties  of  explosions  and  earthquakes,  including 
magnitude,  location,  mechanism,  and  source  time  functions. 
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OBJECTIVE 

The  objective  of  this  study  is  to  develop  and  apply  new  methods  for  the  charac¬ 
terization  of  the  source  of  explosions  and  earthquakes  using  short-period  seismic 
waves,  with  improved  resolution  power  on  source  properties,  including  magnitude, 
location,  mechanism,  and  source  time  functions.  The  accurate  characterization  of  the 
source  with  these  parameters  will  improve  the  regional  monitoring  and  discrimination 
capability  under  a  potential  CTBT  as  well  as  basic  seismological  knowledge  on  source 
processes  involved  in  explosions  and  earthquakes. 


RESEARCH  ACCOMPLISHED 

Our  progress  in  the  past  year  includes  mainly  two  developments:  (1)  including 
Earth’s  aspherical  corrections  in  the  computation  of  synthetic  seismograms  and  in  the 
source  determination  using  recently  developed  three-dimensional  earth  models;  (2) 
developing  a  new  procedure  for  the  suppression  of  noise  and  unmodeled  effects  of 
Earth’s  structure  on  observed  waveforms. 

Although  development  (1)  is  of  importance  to  this  study,  it  is  not  much  different 
from  some  techniques  used  elsewhere  and  thus  is  not  discussed  in  this  report.  We  will 
focus  on  (2),  which  involves  a  new  data  processing  algorithm.  This  algorithm  is 
referred  to  as  the  Earth  simplifying  transformation  (EST),  which  allows  us  to  reduce 
the  complexity  or  to  simplify  observed  waveforms  but  still  retain  the  information 
relevant  to  the  source  determination. 

Since  seismic  waves  recorded  at  regional  and  teleseismic  distances  from  the 
source  are  dominated  by  short-period  surface  waves  in  the  1-100  sec  period  band,  it  is 
desirable  to  use  the  short-period  surface  waves  for  the  determination  of  source  proper¬ 
ties.  However,  this  is  hampered  by  the  complexity  of  the  waves  due  to  the  effects  of 
Earth’s  lateral  and  radial  heterogeneities. 

Earth  simplifying  transformation.  A  procedure,  which  transforms  a  seismic 
signal  into  the  simpliest  form  relative  to  a  given  earth  model.  The  difference  between 
the  simplified  waveform  and  the  original  signal  corresponds  to  the  effects  of  earth’s 
structure  which  are  not  predicted  for  the  earth  model;  and  the  signal  in  the  simplified 
waveform  corresponds  to  the  effects  that  are  predictable  for  the  earth  model  and  infor¬ 
mation  on  the  properties  of  the  source.  This  procedure  enhances  signal-to-noise  ratio 
and  thus  improves  the  resolution  power  of  short-period  surface  waves  in  the  source 
characterization. 

Method.  According  to  Gilbert  and  Dziewonski  [1975]  the  spectrum  of  a  com¬ 
ponent  of  ground  motion  excited  by  a  point  source  at  angular  frequency  co  for  a  given 
earth  model  may  be  given  by 

6 

Uk(x,  CD)  =  £  Ylufx,  xs,  CD)  fj(CO).  (1) 

P=1 

where  uk  is  the  k-th  record  in  a  set  of  seismograms,  with  the  receiver  at  position  x  and 
the  source  at  xs;  are  excitation  kernels  and  f;  represent  six  independent  components 
of  the  moment-rate  tensor.  For  another  event  located  at  x'  with  moment-rate  tensor  ff, 
the  spectrum  may  be  given  by 
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(2) 


Uk(x,  CD)  =  X  ^(X,  x',  CO)  f'(CD). 

Taking  into  account  of  our  imperfect  knowledge  of  the  earth  structure  and  noise 
generated  from  other  sources,  the  spectra  of  observed  ground  motion  for  these  events 
may  be  expressed  as 

Uk(x,  co)  =  uk(x,  co)  ak(x,  co)  +  e(x,  co).  (3a) 

Uk(x,  co)  =  uk(x,  co)  ak(x,  co)  +  e'(x,  co).  (3b) 

where  a  and  a'  represent  effects  of  the  deviation  of  the  earth’s  structure  from  the  earth 
model.  In  the  following  analysis  the  noise  terms  (e  and  eO  are  ignored  and  the  event 
at  xs  is  referred  to  as  the  primary  event  and  the  event  at  x'  as  the  secondary  event. 

In  the  conventional  algorithms  the  moment-rate  tensor  of  the  primary  event,  fj,  is 
determined  from  data  for  the  event,  Uk,  by  solving  the  following  equation 

Uk(x,  CD)  =  X  Yki(x>  XS’  G>)  (4) 

i=l 

which  is  correct  to  the  zeroth-order  in  terms  of  a,  since  implicit  in  the  algorithms  is 
the  assumption 

ak(x,  CD)  =  1.  (5) 

The  question  that  we  shall  address  in  this  study  is  as  follows.  Given  a  set  of 
seismograms  from  two  events,  Uk  and  Uk,  which  are  located  at  a  close  distance  along 
with  the  synthetics  for  the  secondary  event,  uk,  predicted  for  a  given  earth  model,  is  it 
possible  to  determine  the  moment-rate  tensor  of  the  primary  event,  fj?  For  this  pur¬ 
pose  we  define  the  following  transformation 

Up  (X,  CD)  =  Uk(x,  CD)  Uk(X,  CD)  /  Uk(x,  CD).  (6) 

which  is  referred  to  as  the  earth  simplifying  transformation  (EST)  in  the  following 
analysis. 

In  general,  a  and  a'  in  (3)  are  slowly  varying  functions  of  the  source  location 
and  the  moment-rate  tensor,  with  the  principle  term  of  their  Taylor  expansion  being 
independent  of  the  source  location  and  moment-rate  tensor  and  the  first  and  higher 
order  terms  much  smaller  than  the  principle  term.  Therefore  we  assume 

ak(x,  cd)  =  ak(x,  co).  (7) 

Then  the  right  hand  side  of  (6)  becomes  uk(x,  co),  thus  the  spectrum  of  the  EST 
seismogram  may  be  expressed  as 

Up(X,  cd)  =  X  Yki(x>  xs>  ©)  fi(to)-  (g) 

i=l 

which  is  correct  to  the  first-order  in  terms  of  a  and  will  be  used  to  determine  the 
moment-rate  tensor  (f;).  For  the  simplicity  of  our  analysis  we  adopt  a  frequently  used 
assumption:  fj  are  considered  to  be  independent  of  frequency  except  for  a  correction 
for  an  assumed  duration  of  the  source  and  are  regarded  as  the  moment  tensor. 
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Tests.  Several  tests  were  made  using  data  from  earthquakes  located  in  southern 
California  recorded  at  Global  Seismographic  Network  (GSN)  stations.  Figure  1  shows 
the  location  of  the  primary  event  (event  I:  the  January  17,  1994  Northridge  mainshock. 
Mg  6.8)  and  secondary  events  (event  II:  an  aftershock  of  the  1994  Northridge  earth¬ 
quake  sequence,  Ms  6.0;  and  event  HI:  the  June  28,  1991  Sierra  Madre  earthquake,  Ms 
5.2)  along  with  source  mechanisms  for  these  events  obtained  previously  by  various 
investigators.  For  the  primary  event,  three  solutions  of  the  source  mechanism  that  are 
shown  in  Figure  1  are  examined  using  EST  seismograms. 

Figure  2  shows  the  waveforms  of  ground  motion  recorded  at  station  WMQ  and 
synthetics,  bandpass  filtered  between  7  and  30  mHz.  There  is  considerable 
discrepancy  between  the  observed  waveforms  (solid  lines  under  the  labels  2k,  Ri,  and 
Ti)  and  the  synthetics  (dotted  lines  under  the  labels  Zira,  R"1,  and  T")  for  the  primary 
event.  In  clear  contrast,  the  fit  of  the  synthetics  to  the  EST  waveforms  (solid  lines 
under  the  labels  Z®,  Rim,  and  Tin)  is  excellent 

Figure  3  shows  the  transverse  components  of  the  EST  seismograms  along  with 
synthetics  for  the  primary  event  computed  for  the  three  source  mechanisms  that  are 
shown  in  Figure  1  for  several  stations  in  the  azimuths  between  3°  and  12°.  For  each 
station,  synthetics  for  different  source  mechanisms  show  larger  differences  in  ampli¬ 
tudes.  The  amplitudes  of  the  synthetics  for  the  mechanism  from  the  regional 
waveform  inversions  (B  in  Figure  3)  are  in  the  best  agreement  with  the  amplitudes  of 
EST  seismograms  in  comparison  with  other  two  mechanisms  (A  and  C  in  Figure  3). 

Figure  4  shows  the  vertical  components  of  the  EST  seismograms  along  with  syn¬ 
thetics  computed  for  these  source  mechanisms  for  stations  at  various  azimuths.  The 
data  shown  here  represent  only  a  small  portion  of  the  whole  data  set;  however,  they 
are  more  sensitive  to  the  source  mechanisms  discussed  here  and  have  higher  signal-to- 
noise  ratio  than  data  for  other  stations.  Analysis  of  these  EST  seismograms  yields  the 
same  conclusion  as  our  analysis  of  the  transverse  components:  in  general,  the  ampli¬ 
tudes  of  the  synthetics  for  the  mechanism  from  the  regional  waveform  inversions  are 
in  the  best  agreement  with  the  amplitudes  of  EST  seismograms  in  comparison  with 
other  two  mechanisms. 


CONCLUSIONS  AND  RECOMMENDATIONS: 

The  EST  seismograms  provide  superior  resolution  power  in  the  source  characteri¬ 
zation  than  original  seismograms.  This  method  works  better  for  stations  at  close  dis¬ 
tances  from  the  source.  It  is  desirable  to  apply  this  procedure  to  data  from  several 
underground  nuclear  explosions  in  Lop  Nor,  China,  to  characterize  the  source  proper¬ 
ties  of  the  May  21,  1992  explosion.  Another  application  is  to  examine  the  possibility 
to  use  earthquakes  nearby  as  the  secondary  events  to  identify  small  explosions.  The 
application  may  start  with  the  examination  of  the  seismicity  near  several  nuclear  test 
sites  in  the  Central  Asia,  with  both  the  earthquake  and  explosion  being  recorded  at 
regional  and/or  teleseismic  distances. 
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Fig.  1.  Location  (plusses)  and  mechanisms  of  the  1994  Northridge  (I:  mainshock,  II: 
aftershock)  and  1991  Sierra  Madre  (III)  earthquakes.  The  mechanism  for  event  II  is 
shown  for  the  solution  from  the  regional  wave  inversions  and  for  event  ID  for  the 
solution  from  P-wave  first  motions.  Fist  motion  data  are  plotted  with  the  uncertainty 
given  by  the  size  of  each  symbol. 
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Fig.  2.  Data  and  synthetics  for  station  WMQ  (Urumqi,  China).  Symbols  Z,  R,  and  T 
with  subscripts  I,  II,  and  in  indicate  vertical  (Z),  radial  (R),  and  transverse  (T)  com¬ 
ponents  of  ground  displacements,  respectively,  from  various  events  (subscripts  I,  n, 
and  HI  identify  the  events).  Symbols  Z,m  and  R,ra  indicate  the  EST  seismograms 
obtained  using  event  in  as  the  secondary  event,  with  the  synthetics  for  the  event  com¬ 
puted  using  the  focal  mechanism  obtained  from  P-wave  first  motions.  Symbol  T,“  indi¬ 
cates  the  EST  seismograms  obtained  using  event  H  as  the  secondary  event,  with  the 
synthetics  for  the  event  computed  using  the  focal  mechanism  from  regional  waveform 
inversions.  Dotted  lines  indicate  synthetics  for  the  primary  event  computed  using  the 
focal  mechanism  from  P-wave  first  motions.  Amplitude  is  arbitraryly  set  to  show 
waveform  coherence.  Time  is  measured  from  the  origin  time  of  the  corresponding 
event. 
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Fig.  3.  Transverse  components  of  EST  seismograms  (T,n)  at  various  stations  obtained 
with  the  synthetics  for  event  II  computed  using  the  focal  mechanism  from  the  regional 
wave  inversions,  with  amplitudes  being  given  in  parentheses  in  units  of  microns 
(amplitudes  of  T,"  with  the  synthetics  computed  for  event  II  using  the  focal  mechanism 
from  P-wave  first  motions  are  shown  in  brackets).  Symbols  A,  B,  and  C  indicate  syn¬ 
thetics  computed  for  the  primary  event  using  a  seismic  moment  of  1.2  x  1019Nm  and 
various  focal  mechanisms:  first-motion  solution  (A),  regional-wave  solution  (B),  and 
CMT  solution  (C). 
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Fig.  4.  Vertical  components  of  EST  seismograms  for  various  stations.  See  caption  to 
Figure  3  for  symbol  idetification. 


701 


